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We calculate the electronic and optical excitations of polythiophene using the GW approximation 
for the electronic self-energy, and include excitonic effects by solving the electron-hole Bethe-Salpeter 
equation. Two different situations are studied: excitations on isolated chains and excitations on 
O , chains in crystalline polythiophene. The dielectric tensor for the crystalline situation is obtained 

by modeling the polymer chains as polarizable line objects, with a long-wavelength polarizability 
tensor obtained from the ab-initio polarizability function of the isolated chain. With this model 
dielectric tensor we construct a screened interaction for the crystalline case, including both intra- 
and interchain screening. In the crystalline situation both the quasi-particle band gap and the 
exciton binding energies are drastically reduced in comparison with the isolated chain. Ifowever, 
. the optical gap is hardly affected. We expect this result to be relevant for conjugated polymers in 

general. 

^ '. 78.40.Me,7I.20.Rv,42.70.Jk,36.20Kd 
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' I. INTRODUCTION 

^ ' Semiconducting conjugated organic polymers have received increasing interest in recent years, especially since the 
, discovery of electroluminescencelj of these materials. The charge carriers and excitations in these materials have 
■ been studied extensively both experimentally and theoretically, but many important fundamental issues still remain 
, unsolved. For instance, the magnitude of the exciton binding energy in these materials is still disputed.a This is a 
very important quantity, since e.g. in photovoltaic devices (solar cells) one would like to have a small binding energy, 
which facilitates the fast separation of charges, while in electroluminescent devices such as LEDs a larger exciton 
. binding energy, to increase the probability of fast (radiative) annihilation of electron-hole pairs, is desirable. 

In conventional semiconductors such as Si and GaAs the optical excitations are well described in terms of very 
weakly bound electron-hole pairs (so-called Wannier excitons) with a binding energy of the order of 0.01 eV. In 
S ' crystals made of small organic molecules such as anthracene, the exciton is essentially confined to a single molecule 
(Frenkel exciton), leading to a binding energy of the order of 1 eV. The question is where exactly conjugated polymers 
(— ( fit in bet"WKen conventional semiconductors on the one handpaiad molecular crystals on the other: negligibly small (0.1 
Q ' eV or lessQ), intermediate (~ 0.5 eVcl), and large (^1.0 eVQiD) binding energies have been proposed. 
O ' Ab-initio calculations, on a variety of conjugated polymers, within the Local Density Approximatipa of Density 
^ [ Functional Theory (DFT-LDA) yield equilibrium structures in very good agreement with experiment Bt12I The Kohn- 
Sham gaps in these calculations are typically 40% smaller than the optical band gap (absorption gap). In cases where 
calculations for both isolated chains and the crystalline situa±ixm were performed, small differences (0.1 to 0.2 eV 
for gaps of ~ 3.0 eV) in Kohn-Sham band gaps were found.EjllJ However, it is well known that the Kohn-Sham 
eigenvalues formally cannot be interpreted as excited state energies.Ej Moreover, excitonic effects are not taken into 
account in these calculations. j— , 

An ab-initio many-bpdy calculation within the GW ApproximatioiiE3 (GWA) was performed for poly-acetylene 
(PA) by Ethridge et aZ.EJ They claim that their quasi-particle (QP) gap, excluding excitonic effects, is in agreement 
with the epieerimental absorption gap. This result seems to be in contrast with a more recent calculation by Rohlfing 
and Louial3 of both one- (QP) and two-particle (exciton) excitation energies for PA and poly-phenylene-vinylene 
(PPV) chains. Their absorption gaps are in good agreement with experiments, but the inclusion of excitonic effects 
proves to be crucial for this. However, their exciton binding energy of 0.9 eV for PPV is much larger than recently 
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obtained experimental values: 0.35 ±0.15 eV for an alkoxy-substituted PPV,i and 0.48 ±0.14 eV0 for unsubstituted 
PPV. rn 

In a recent Letter,tZI hereafter referred to as I, we focused on the differences in excitations between an isolated 
polythiophene (PT) chain, see Fig. ^, and crystalline polythiophene. For the isolated chain, we found an absorption 
gap in good agreement with experiment, but the energy differences between the various exciton levels were too large.li^ 
After including the screening by the surrounding chains, both the optical gap and the exciton transition energies were 
in good agreement with the experimental values. The difference in screening between an isolated polymer chain and 
a condensed polymer medium can be explained as follows. For jaii isolated quasi-one dimensional system, such as a 
single polymer chain in vacuum, there is no long-range screening.tj lII A way to understand this is to realize that if we 
have two charges on a polymer chain with a separation larger than the width of the chain (~ 7 a.u.), most field lines 
connecting the charges will be outside the chain. If the chain is embedded in a medium (possibly, but not necessarily, 
consisting of similar chains), the medium will provide a long-range screening of the Coulomb interaction. The screened 
Coulomb interaction determines both the QP energies (including the band gap) and the exciton energies. Long-range 
screening reduces both the QP band gap and the exciton binding energies. Apparently, there is a near cancellation 
between the change of the QP gap and the exciton binding energy, meaning that the optical absorption gap, which is 
the difference between the two, is influenced much less by introducing screening. 

In I, a method for the calculation of the dielectric tensor of crystalline polythiophene from the ab-initio single chain 
polarizability function was introduced, without giving any details. Further, some novel technical procedures in the 
GWA calculation were used, in particular in the handling of the Coulomb divergence both in real and reciprocal space. 
Details of these approaches, as well as of the calculation of the quasi-particle energies and exciton binding energies, 
will be explained here. The paper is organized as follows: in the next section (||) we explain the computational 
methods employed to calculate the quasi-particle bandstructure, to regularize Coulomb interaction, to calculate the 
exciton binding energies and the dielectric tensor. In Section [II, we will present results for the electronic and optical 
excitations of both the isolated chain and bulk PT. In section [V we will discuss these results and compare them to 
other calculations, and draw our conclusions. 



II. COMPUTATIONAL METHODS 



Many successful ab-initio calculations of the QP band structure of conventional anorganic semiconductors have 
been performed within the GW^AI13 for the electronic self-energy E of the onepQaxticle Green function. Very recently, 
progress has been made in the evaluation of the two-particle Green fuiiction,E3'Ej from which the optical properties 
can be obtained. This is done by solving the Bethe-Salpeter equationE3'E3 (BSE), which can be mapped onto a two- 
body Schrodinger equation for an electron and a hole forming an exciton. We will use these approaches to calculate 
the QP band structure and exciton binding energies of PT. The calculational scheme is as follows: first we perform a 
DFT-LDA-bascd Car-Parrinello calculation, from which we obtain atomic positions, wave functions and ground-state 
energies. We use these as input for the GWA calculation, which yields the QP excitation energies. With the DFT- 
LDA wave functions and the QP energies we calculate the two-particle excitations by solving the BSE. This scheme 
is first applied to the single chain and next to the crystal. We assume the same atomic geometry for the ground and, 
excited states, i.e. the coupling between electronic and lattice degrees of freedom is neglected. ExperimentaJ datatS 
indicate that energy shifts due to lattice relaxations are of the order of 0.1 eV in PT. DFT-LDA calculation^ predict 
a hole-polaron relaxation energy of 0.04 eV for 16T (nT is an oligomer consisting of n thiophene- rings). A similar 
calculation predicts a triplet exciton relaxation energy of 0.2 cV for 12T.E3 Singlet relaxation energies are typically 
smaller. These values, calculated for oligomers, are upper bounds for the values in the polymer, since in the oligomers 
the excitation is confined, leading to a larger local deviation from the ground state density and hence to a larger 
relaxation energy than in the polymer. 



A. The quasi-particle equation 

We start our calculations with a pseudo-potential plane-wave DFT-LDA calculationi of a geometry-relaxed PT 
chain in a tetragonal supercell. The plane wave cut-off energy is 40 Ry. The length inutile chain direction ax was 
optimized and found to be — 14.80 a.u. (experimental values range from 14.65 a.u.E3 to 15.18 a.uB3). In the 
perpendicular directions we found that a separation of ay = a^ = 15.0 a.u. is enough to consider the chains in the 
DFT-LDA calculation as non-interacting. The two rings in the unit cell are found to be co-planar and we choose 
them in the y = z plane. We use Hartree atomic units (with the Bohr radius ag as unit of length and the Hartree as 
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unit of energy) throughout this article, unless specified otherwise. The one-particle excitation energies are evaluated 
by solving the QP equation: 
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Mr) + / [Vppir,r') + S(r,r',S„fe)] </)„fc(r')dV - EnkM^), (1) 



where Vh is the Hartree potential, Vpp the non-local pseudo-potential of the atomic core, and S the electronic self- 
energy. Since in practice the DFT-LDA wave functions and the QP wave functions are almost identical, we use the 
former in all calculations. In DFT-LDA, S is approximated by: 

^{r,r',Lo) = V^c{r)S{r-r'), (2) 

where Vxc is the exchange-correlation potential of the homogeneous electron gas. In the GWA, S is approximated by 
the first term of the many-body expansion in terms of the one-particle Green function G and the screened Coulomb 
interaction W of the_Kystem.E^ In order to calculate IL we follow the real-space imaginary-time formulation of the 
GWA of Rojas et al.jEi in its mixed-space formulation.L3 In this formulation, we transform non-local functions F{r, r') 
to functions Ffc(r,r') 

Ffe(r,r') = ^^^(r-f na,i;,r')e-^'=("+"'^^-^'), (3) 

n=l 

F(r,r') = — 5]Ffc(r,r')e^'=(--^'), (4) 

where Uk is the number of equidistant /c-points in the ID Brillouin zone (BZ). We use periodic boundary conditions: 
F{r + nka^x, r') = F{r, r' -t- Ukaxx) = F(r, r'). The functions Ffe(r, r') are fully periodic: Fk{v + axX^ r') = -Ffe(r, r' + 
axx) = Ffc(r,r'), so that r and r' can be chosen in the unit cell. We calculate the one-particle Green function for 
imaginary times 

^ u„fe(r)<fe(r')e-(^--^-)^ for r < 0, 

Gfc(r,r',zr) = <( ^ (5) 
-i Ucfe(r)u*fe(r')e-(<^<^'--'^)^ for r > 0, 

c 

where ■(i„fc(r) — (/)„fe(r)e^*'^^ and e„fc (with n = c,v) are DFT-LDA wavefunctions and the corresponding energies 
(with V and c referring to valence and conduction states, respectively), ep is the Fermi energy (set in the middle of 
the DFT-LDA gap). We further calculate the irreducible single-chain polarizability function in the Random Phase 
Approximation (RPA) 

Pk{r,r',iT) = -2z^G,(r,r',zT)G,_fc(r',r,-iT), (6) 

9 

the screened Coulomb interaction 

Wk{r,r',iLu)= [i-\r,r')~Pk{r,r\iLu)y\ (7) 

where Vk(j,r') is a cut-off Coulomb interaction in mixed space, discussed in the next Section, and we calculate the 
electronic self-energy 

I]fe(r, r', it)=iJ2 Gg(r, r', iT)Wk^,{r, r', zr). (8) 

We calculate all the above two-point functions on a double 24 x 24 x 24 real-space grid for r and r' in the unit cell. This 
corresponds to a plane wave cut-off of 25 Ry. The total number of valence and conduction bands taken into account 
was 300. In Eqs. (^), (0) and (^) we switch between time- and frequency domain using Fourier transforms. Our 
imaginary-time grid has an exponential spacing (0.25 a.u. near r = 0, up to a spacing 6 a.u. near Tmax — 32.0 a.u.) 
and we interpolate to a linear grid when using the Fast Fourier Transform (FFT) to imaginary frequency. A similar 
exponential grid is used for the imaginary frequency. We split the self-energy in an exchange part and a correlation 
part S=: 
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E^(r, r') = ~^S)v,^k{r', r), (9) 

S^(r,r',zT) =^^G,(r,r',^r)W^-'■,(r,r,^T), (10) 

where (5 is an infinitesimally small positive time, and W^'^'^ is the screening interaction: 

W^"ir,r',iu;) = Wk{r,r',iuj) - Vk{r,r'). (11) 

In the calculation of we use a ID Brillouin-zone sampling of 10 equally spaced fc-points, and in the calculation 
of S'^ 4 fc-points (since the screening interaction is short-ranged, the convergence of with respect to the number 
of fc-points is faster than that of S^). With the above parameters, the calculated QP gap has converged to within 
about 0.05 eV. From a two-pole firt-pn the imaginary-frequency axis an analytical continuation to the real-frequency 
axis is obtained: E'^(ia;) E'^(w)B2l Subsequently, the QP equation Eq. (|l|) can be solved by replacing the QP wave 
functions by the DFT-LDA wave functions and obtaining Enk iteratively: 

Enk = + {(bnk l^UEnk) + ~ V^'l 4>nk) ■ (12) 



B. Treatment of the Coulomb interaction 



We have developed a novel procedure to deal with the k = and r = r' singularities of the Coulomb interaction 
?;k(r, r'). We will first describe the procedure for a 3D system and later explain the specific adaptations of this 
procedure we used for our quasi- ID system. 

In reciprocal space the Coulomb interaction is given by: 

Att 

«k(K,K') = ^^^<5K,K'. (13) 

We replace fo(0,0), which would be infinite in Eq. ( |l3|) by a finite value, which is obtained in the following way. We 
evaluate the integral over all space of the Coulomb interaction multiplied by a Gaussian: 

/a= / d^q^e-"'^ = 8^^/^, (14) 
J \J a 

and evaluate the corresponding sum, excluding the singularity for K = k = 0: 

= AV 5] V(K, K)e-"l'^+K|' , (15) 

k,K 

where k g IBZ, the first Brillouin zone of the 3D lattice. AV is the volume per (k, K)-point and the prime indicates 
that k = K = is excluded in this sum. We now put: 

Vk=oi'K = 0, K' = 0) = lim [/„ - S^] . (16) 

Finally, we obtain v{r — r') by a discrete FFT of v^CK, K') to real space. We find a finite value for i;(r — r' = 0), 
solving at the sameJime the problem with the Coulomb singularity for r — r' = 0. In the original formulation of the 
space-time method,Lil the authors used a grid for r' offset with respect the r-grid in order to avoid this singularity. 

In order to study a truly isolated chain, which is a quasi-lD system, we have to avoid 'crosstalk' between periodic 
images of the chain in the perpendicular directions. We do this by dividing space into regions of points that are closer 
to the atoms of a specific chain than to those of any other. Subsequently, we cut off the Coulomb interaction u(r — r'), 
obtained in the way described above, by setting it zero if r and r' belong to different regions. Thus we obtain an 
interaction ^(r, r'). In the construction of the Coulomb interaction v(r — r'), we take a regular grid of k-points with a 
spacing in the y- and z-direction approximately equal to that in the x-direction. From the cut-off interaction ^(r, r') 
we obtain Vk{v,r') in mixed space from Eq. (^ with k now in the ID Brillouin zone. 
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C. The Bethe-Salpeter equation 



The two-body electron-hole Schrodinger equation related to the BSE is solved by expanding the exciton wave 



functions <i>(re,r;i) in products of conduction (pcki'^e) and valence wave functions 4'vki'^h) 
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$(re,r/,) = ^ Afcc„</>cfc(re)0*fe(r,,). (17) 

k,c,v 

Here we have restricted our discussion to excitons which have zero total momentum, since only these are optically 
active. As we are interested in the lowest lying excitons, an expansion in the highest occupied valence (tt) and lowest 
unoccupied conduction (tt*) bands is sufficient to converge the exciton energies to within 0.1 eV; energy differences 
are converged even better. Below, we will give all energies in aath a precision of two decimal places. The exciton 
binding energies E}, follow from the Schrodinger-like equationiEacil 

[Eck — Eyk — Eg + Ef,]Akcv + ^ ['^^kcv,k'c'v'^s,0 — Wkcv,k'c'v']Ak'c'v' = 0, (18) 

k' c'v' 

where Eg is the QP band gap. Eh the exciton binding energy and Wkcv.k' c'v' are the matrix elements of the static 
{lj ~ 0) screened interaction 

Wkcv^k'c'v ^ Jd'^rJ d3rV„fc(r')Cfc(i-)W^(r,r',L^ = 0)C'fc' (r')0c'fc' (r), (19) 
and V^^^ k'c'v' exchange matrix elements (present for singlet excitons, s = 0, only) of the bare Coulomb interaction: 

VScv,k'c'v' = ld?v [ d3r>;,,,(r)0:,(r')«(r,r')0.fc(r')</>c'fc'(r). (20) 



The integrals over space in Eqs. (y^ and (|20|) are in the calculations replaced by summations over our real-space 
grid. We use wave functions and energies on a grid of 8 fc-points and extrapolate to a grid of 100 fc-points. 

Formally, dynamical screening effects may only be ignored in the BSE if Eg ^ E},. However, since it has been 
showrL±hat dynamical effects in the electron-hole screening and in the one-particle Green function largely cancel each 
other ,Ej this approximation is nevertheless valid, even if the relation Eg ^ Ef, does not strictly hold. 

We calculate an approximate exciton size Oex by fitting the exciton coefficients A^cv to the hydrogen-like form: 

(l + a2,fc2)2- ^''> 

Note that in fact the exciton is highly anisotropic. Nevertheless, Eq. ( ^ ) gives pretty good fits and can be used to 
get a qualitative impression of the (relative) size of the excitons. 



D. Inclusion of interchain screening 

As mentioned in the introduction, in a quasi-lD system, such as an isolated chain of a polymer in vacuum, there is 
no long-range screening. For a meaningful comparison of our calculations to the experimental data, which are obtained 
from either films or bulk polymer material, both the intra- and the interchain screening are important, and only the 
latter is long-ranged. It would be desirable to perform a GWK and exciton calculation for a 3D crystal structure 
of PT, but the amount of computational work involved is as yet prohibitively large. Since PT samples prepared in 
many different ways show very similar optical behavior, we expect the details of the interchain screening not to be 
extremely important. This consideration leads us to propose the following approximation for the interchain screening 
interaction, defined analogously to Eq. (pj]): 

W,Xr(r,H = (1 - e-/'---) { [el{^u;)x^ + e\\{^u^)e ^{zLo){y^ + z^)Y"^ - v{v)] , (22) 

where e_\^{iuj) and e\\{iu}) are the ab-initio frequency-dependent dielectric constants perpendicular and parallel to the 
chain, respectively. The counter-intuitive combination of dielectric constants and coordinates in Eq. (^2|) results 
from solviMjthe Laplace equation for a point charge in a homogeneous, anisotropic medium with dielectric constants 
£|| and e_LHj The prefactor takes care of a smooth cut-off for distances smaller than ri„tcr, for which the interchain 
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screening should be replaced by the intrachain screening. Eq. (g2|) has the correct behaviour for distances larger than 
the interchain distance rintcr, for which we take 10 a.u., which is typical for the experimental crystal structures of 



Refs. E9 and 30 



The total screened interaction for the bulk system then becomes: 

W^totai(r, r', zuj) = W^Z.ir, r' , ilo) + W^^^U^ - r', ioj) + v{r - r'), (23) 

where Wi^jj"j.3^ is the intrachain screening already calculated with Eq. (0). The screened interaction Wtotai is correct 
at short range, where the interchain screening is vanishingly small compared to the 1/r-divergence of the intrachain 
screening, and at long range, where the intrachain screening vanishes due to its quasi-lD nature. Of course, for 
intermediate ranges, it is not strictly allowed to simply add the parts representing long- and short-ranged screening, 
but we expect Eq. (|2|) to give a reasonable interpolation there. 

Note that the interchain screening part given by Eq. ( p^ ) is long-ranged by construction, and 8 fc-points are now 
needed to converge the corresponding self-energy Sf„ter from Eq. (p^. On the other hand, the required number of 
real-space grid points in order to calculate ^f^tcr l^^s than before, because Wj'^tgj. is a very smooth function of r; a 
12 X 12 X 12 real space grid turns out to be sufficient. The total self-energy can be expressed as: 



■^total 



— ^hitra + ^hitcr + (24) 



Because the self-energies in this equation are additive, we can reuse the self-energies Sf^tra which we have 

already calculated for the isolated chain. 

The overlap between wave functions, and therefore the electronic coupling between neighboring chains, is very small. 
This means that we can use the isolated-chain wave functions to calculate the Green function and self-energy. This 
obviously implies that in our exciton calculations we restrict ourselves to excitons in which we take the electron and 
hole are on the same chain (so-called intrachain excitons). In summary, the only, but essential, difference between our 
calculations for the isolated PT chain and bulk PT is in the use of an interchain screened interaction. 



E. Dielectric tensor of crystalline PT 



In order to construct the screened interaction of Eq. (|22), we have to calculate the dielectric tensor of bulk PT. We 
do this for the crystalline structure of Ref. |2^, which is reproduced in Fig. |^. We use a model in which the chains are 
replaced by polarizable line objects with a polarizability tensor obtained from the single-chain polarizability function. 
The principal axes of the chain are the following: 



XI =x, X2 = -^(5 -I- y), X3 = -^(z - y). 
The full polarizability function X(r,r',ia;) of a single chain is given by: 



X{r, r', ioj) = P(r, r', iui) 



dr"dr"'P{r, r", icujWintUr" , t'" ,iuj)P{t"' , r', lu) 



X(°)(r,r',iu;) +X(i)(r,r',it^). 



(25) 

(26) 
(27) 



The long- wavelength (g — s- 0) polarizability tensor x P^r unit of chain length of a single chain in the {xi^X2,xz) 
coordinate system is diagonal and has diagonal elements given by: 



Xi(iw) = lini 



and for j=2,3: 



Xj(iw) 



drdr' XjX{Y. r', iujjx'^. 



(28) 



(29) 



The calculation of Xjii'^) has been performed with 4 fc-points, with the exception of Xi'\ for which it proved to be 
necessary to use 8 fc-points. 

If we now approximate the chains by polariMkble line objects with the above polarizability tensor, we can calculate 
the macroscopic dielectric tensor of the crystalE^ of these chains. This is done by a procedure of which the details are 
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given in the Appendix. The axes of the crystal unit cell are denoted by a, b and c = x = xi. Dropping the frequency 
dependence in the notation, we find the following expression for Sc- 

ec = l + — ^, (30) 
where A is the surface area per chain in the plane perpendicular to the chain. For and we find: 

where j — a,b and x-y is the effective polarizability of the chain along the 7-axis. In the Appendix details of the 
calculation of Ea, £b, and Sc are given. 

To retain the tetragonal symmetry in our calculation (in order keep the computations feasible), we average ea{ii^) 
and ef,(i(jj), which are not very different, to obtain £j_(ia;). For e\\(iuj) we take £c(*t^)- Note th at fo r using the screened 
interaction of Eq. ( p2| ) to in the implementation of the GWA formalism presented in Section p^I A| , we have calculated 
the dielectric constants ej_{iLu) and £||(ia;) along the imaginary frequency axis. 



III. RESULTS 



A. Isolated chain 

The calculated GWA QP band structure (together with the DFT-LDA band structure) is shown in Fig. ^ left 
panel. We find a minimal band gap Eg at P of 3.59 eV, which is quite large compared to the DFT-LDA value of 1.22 
eV. The effective masses, m* = l/h^{d^E/dk^)~^, of the tt and tt* bands at P, which are 0.15 and 0.17 rUe (with rUe 
the free electron mass) in DFT-LDA, are reduced by about 15% in the GWA to 0.13 and 0.15 me- This corresponds 
to an increase of the band width from 1.96 and 1.51 eV in DFT-LDA to 2.48 and 1.81 eV in the GWA, for the tt and 
TT* bands, r*espectively. In an earlier GWA study, a similar increase of the bandwidth was found for a wide variety of 
materials .£3 

The lowest-lying singlet exciton (^Bu) has a binding energy Ei, of 1.85 eV. The size Oox of this exciton, calculated 



using Eq. (|21 
^'(re,r/i), we 
the electron, 



), is 12 a.u., i.e. less than two thiophcne rings. To give an impression of the exciton wave function 
lave plotted in Fig. (top panel) the probability to find the hole at a distance x/j along the chain from 



Prob(a:,,)^ dyhdzhMre,rh)\', (32) 



where the electron coordinate re is taken 1 a.u. from the inversion center, in the direction perpendicular to the polymer 
plane (for the electron coordinate in the inversion center, this probability would be zero due to symmetry). We have 
plotted Prob(a:?i) for both the ^Bu and ^Ag excitons. 

As thejaptical gap is given by Eo = Eg — Ei,, we have Eq = 1.74 eV, in good agreement with the experimental value 
of 1.8 eVll3 (see Table ||). While there is good agreement for the optical gap, the differenee between the ^Bu and ^Ag 
binding energies of the isolated PT chain is definitely not in agreement with experiment,E3 see Table |. Moreover, the 
^Bu exciton binding energy oLl.85 eV is very large compared to values currently discussed in the literature, which 
range from ~ 0.1 to ~ 1.0 eV.I 



B. Dielectric properties 

We calculate the polarizabilities per unit length Xji^^) with Eqs. ( |2^ ) and (p9|). The obtained w = values are 
listed in Table Note that the polarizability along the chain, i.e. in the direction of the extended carbon 7r-system, 
is much larger than those in the perpendicular directions. This difference is reflected in the dielectric constants 
e~f{iuj) calculated using Eqs. (^0|) and (^T|); the dielectric constant along the chain is much larger than those in the 
perpendicular directions. In real systems with disorder the conjugation length will be finite, which will reduce £||. 
Note, however, that the perpendicular dielectric constant e±_ plays the dominant role in the interchain screening of 
Eq. (p^) along the chain. 
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C. Crystalline polythiophene 



The resulting band structure, calculated using the bulk screening from Eqs. (|2^) and ( p3| ) is given in Fig. right 
panel. The QP gap Eg has decreased to 2.49 eV; the ^B^ exciton binding energy is 0.76 eV (see Table |). Hence, the 
predicted optical ga,p is 1.73 eV, virtually unchanged from the isolated chain results of 1.74 eV and in good iSffreement 
with experiment .tSl Note that the absorption gap of Ref. |l^ is 2.0 eV, also found in earlier work on PT,Ea but the 
luminescence gap is 1.8 eV. There are two reasons why we should compare our result to the latter gap. The first 
reason is that absorption occurs everywhere in a sample, both in the ordered and disordered parts, but luminescence 
occurs predominantly in the most ordered parts with the longest conjugation lengths. Thisj-is because, prior to 
recombination, excitons diffuse to those parts of the sample where they have the lowest energy.EZl The second reason 
is that after photoexcitation, the rings, which may be twisted around their common C-C bond, tend to co-planarize 
in the excited state, due to the fact that the excited state is slightly more quinoid than the aromatic ground state.El 
As we are performing our calculations for a perfect, co-planar chain of FT, we should therefore compare our optical 
gap to the luminescence gap. Note that in principle it is possible that excitons trapped in defects or disordered parts 
of the sample to have a lower energy than in a fully conjugated, defect-free polymer. However, the luminescence 
spectrum of Ref. |l8| can be fully understood in terms of the exciton decay and its vibronic side bands, which 
means that such defects are either rare or that excitons trapped by such defects decay non-radiatively. 

What is very important, is that the relative exciton energies (also listed in Table ||) are now also in good agreement 
with experiment. The sizes of the excitons have increased by ^ 50%; the size Oex is now 18 a.u., or slightly more 
than two rings. In Fig. [s] (bottom panel) it is clearly seen that the excitons are larger than the corresponding excitons 
on the isolated chain (top panel). 

In order to test the sensitivity of our results to the precise value of the cutoff distance rintor in Eq. (|2^ ) we performed 
similar calculations for rintcr = 8 a.u. and rintcr = 12 a.u. These data are also listed in Table |. The QP band gaps 
are 2.32 and 2.61 eV, respectively. The ^B^ binding energies are 0.64 and 0.86 eV and hence the optical gaps are 1.68 
and 1.73 eV, respectively. This means that the optical gap is quite insensitive to the choice of Hntor- This is consistent 
with the fact that in the limit rintor ~^ oo, which corresponds to no interchain screening, we should find the isolated 
chain absorption gap of 1.74 eV. The energy differences between the excitons are even less sensitive to rinter- The 
good agreement with experiment and the fact that especially the optical gap and the energy separation between the 
excitons are hardly influenced by varying rintcr are also a posteriori justifications for our model screening interaction 
Eq. (H). 



IV. CONCLUSIONS AND DISCUSSION 



Summing up, we have calculated the quasi-particle band structure and lowest-lying exciton binding energies of an 
isolated polythiophene chain and crystalline polythiophene. For the isolated chain (where there is only intrachain 
screening) we find a large band gap and large exciton binding energies, due to the absence of long-range screening. 
After including interchain screening, which is responsible for the long-range screening in bulk polythiophene, we find 
that both the band gap and exciton binding energies are drastically reduced. However, the optical gap is hardly 
affected. We suggest that these conclusions hold for conjugated polymers is. general. 

This sheds light on the fact that the calculations by Rohlfing and Louiella on isolated chains of PA and PPV yield 
good results for the optical gaps, whereas their lowest-lying singlet exciton binding energy of 0.9 eV for PPV is in excess 
of recent experimental values of 0.35 ± 0.15 eV, 3-pbtained by a direct STM measurement for an alkoxy-substituted 
PPV, and 0.48 ± 0.14 eV for unsubstituted PPV.E3 The inclusion of interchain screening effects will drastically reduce 
their calculated binding energy and may well lead to agreement with this experiment. Clearly, it would also be very 
interesting to repeat the experiment in Ref. ^ for polythiophene and polyacetylene. Interestingly enough, a value of 
0.4 eV is obtained for the exciton binding energy in PPV by mepis of an effective-mass appromixation in which the 
electron-hole interaction is derived from a bulk dielectric tensor.ES The difference of about a factor of two in exciton 
binding energy between crystalline PT and PPV can, at least qualitatively, be explained by the differences is, reduced 
masses ^ (l//x = 1/m^ -I- I/tott*) of PT and PPV, for which we find /^^"^ = 0.08me, while /i^^^ = 0.04TOe,L3 both in 
DFT-LDA, and by the fact that in an effective-mass approximation the binding energy is proportional to /i. Of course, 
these arguments, which are qualitative only, do not take away the need for ab-initio calculations on the crystalline 
phase of PPV. |— , , 

Further, the apparent discrepancy of the results for PA by Ethridge et alSj and those of Rohlfing and Louie,t3 
can be understood. The latter find, for an isolated chain, a QP gap of 2.1 eV and an exciton binding energy of 0.4 
eV, yielding an absorption gap of 1.7 eV. The former find a QP gap of 1.86 eV and do not include excitonic effects. 
This calculation, however, is performed for one PA chain in the same volume as a PA chain in a crystal would have. 
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Therefore, this calculation is in fact one for a bulk situation, which means that this QP gap is by our arguments 
expected to be smaller than that of Rohlfing and Louie. Furthermore, our arguments predict an exciton binding 
energy in bulk PA considerably smaller than the 0.4 eV of Rohlfing and Louie. 

We conclude that a correct many-body description of the electronic and optical properties of bulk polymer systems 
should include the effect of interchain screening. An important consequence of this conclusion is that neither Hartree- 
Fock nor DFT-LDA calculations should be relied upon in this context, since Hartree-Fock does not contain screening 
effects at all and since the exchange-correlation potential in DFT-LDA only depends on the local density and cannot 
describe the non-local effects due to the long-range screening. Moreover, since exciton effects play such a large role 
in conjugated polymers, it is essential to take these effects into account. 
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APPENDIX A: CALCULATION OF THE CRYSTAL DIELECTRIC TENSOR WITHIN A LINE-DIPOLE 

MODEL 



We apply an electric field Eappi(r) = Eqc**' '" (and we will take the limit k 0), where Eq and k are parallel to 
the a,b or c-axis of the crystal (sec Fig. |^) to calculate Ea, £& and £c, respectively. The applied field Eappi leads to an 
induced field Ein(j(r); the total microscopic field Eniici(r) is then given by: 

Emicr(r) = Eappl(r) + Eind(r). (Al) 

We define p = ud + vb with ^ + v"^. Note that there are two different chains: the A type, at the corners of the 
unit cell, and the B type at the center of the unit cell. For the A and B chain we have: 

Pa{^)^X_^-KucA^,P^O) (A2) 
Pb{^) ^Xis-'^LcAx,P^^a+h) (A3) 

with Pj^{x) (pg(x)) the long-wavelength dipole moment per unit length of the A (B) chain and (Xj?) polariz- 
ability tensor of the A (B) chain calculated with Eqs. (^8|) and (^) and using the relations: 

X^ = Ua'-X-Ua^ X^^U^'-X-ILb. (A4) 

with C/^ and C/g are the rotation matrices relating the (2:2 ,2:3) coordinate system to the (a, 5) coordinate system 
(c = xi): 

C/^=f'°1f ~''1f-"!V (A5) 

/ cos{^ - a) sm{^ - a) \ , . 

" I - sin(f - a) cos(i^ -a))- ^^^^ 



The prime in Eqs. ( |A2| ) and (AS) indicates that the field caused by the chain itself is excluded. We will refer to our 
model, in which a PT chain is represented by an homogeneous line with a certain dipole moment per unit length p, 
as a 'line-dipolc'. 

In CGS units the dielectric tensor e is defined as: 

E(r) -K47rP(r) = e-E(r). (A7) 

where E(r) is the macroscopic field, and P(r) is the macroscopic polarization. For each direction of the applied field, 
we will calculate Eind(r), evaluate the macroscopic fields E(r) and P(r) by averaging, and solve Eq. (A7) to obtain 
the dielectric tensor e. 
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1. Calculation of £c 



For Eappi and k parallel to x (and hence to c and also xi), we have for both the A and B chain from Eqs. ( [A^ ) 
and ([A3|) : 



The field induced by a line-dipole on the a;-axis is given by: 

Pxe'-^'^' {x ~x') , 



Eind(r) = -V$(r) = -V 



r — r' 



'13 



(A8) 



(A9) 



where $ is the electrostatic potential and we have used the fact that Pxi^') — PxS ■ Evaluation of Eq. (A9) yields: 



Si„d,.(r) = -2fcV^o(pfc)e^'=^ 



(AlO) 



where Kq is a zeroth order Bessel function of the third kind. From here on, we omit the factor e'*^^. We can calculate 
the total microscopic field at the a;-axis, due to both applied and induced fields, for a crystal of line-dipoles, by 
summing over all line-dipoles but the one at the origin: 



(All) 



where the positions of the other chains are given by p^. In the limit fc — > 0, we can replace the sum by an integral 

2Trp, 



lim^i?i„d,.(pl) = ^ / 



p'dp'Koip') 



47r 



(A12) 



(A13) 



where p' = pk an d A — ab/2 is the area of the two dimensional unit cell per chain. Substitution of Eq. ( [A13| ) and 
( [aTiI ) in Eq. yields: 



A + 47rxi ' 



-E: 



appl,x- 



Since = Px/ A^ we have: 



P. 



A + Anxi 



appl,a; ■ 



(AM) 



(A15) 



The m acroscopic field Ex is the average over the two-dimensional unit cell of the microscopic field as given by Eq. 
(All) for general p, but now including the chain at = in the sum: 



Ex = i^appi,x -I- lim -5- / py^EindAP ~ Pi) 

fe^O ZA coll ^ 

= ^^appi,x -I- lim — / p^pV" Si„d,x(p - Pi) 

k^O A J '-^ 



E: 



appLx 



A 



Combining this with Eqs. (A7) and (A15) we obtain Eq. 



Sr. = 1 



47rxi 
A 



(A16) 
(A17) 
(A18) 

(A19) 
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2. Calculation of and st 



We now take Eappi(r) and k parallel to a. The derivation for Eappi(r) and k parallel to b is equivalent. The dipole 
moments of the chains must satisfy Eqs. ( |A2| ) and (|A3|). The field induced by the chain at the origin is given by: 



Eind(r) = -V$(r) - -Vy" 



where 



r — r 
Mp) ■ Pa 



(iu-' 2 4uv 



,/|3 



dx' 



(A20) 
(A21) 

(A22) 



in the two dimensional (a, 6) coordinate system (the dipole moment in the c-direction is zero and hence we work with 
2x2 instead of 3 x 3 matrices) . The microscopic electric field £1^^;^,^. at the origin, excluding the field induced by chain 
at the origin itself, is given by: 



where 



E;„icr(P = 0) = Eappi(p = 0) + • +iy:6 • Pb, 



M _^ = lim M{pi) cos kui, 

~* PieA,pi=io 



M p = hm AL{Pj ) cos ki 



These sums are evaluated in the next subsection. Substitution of Eq. ( A23| ) in Eq. (A2) and solving yields: 



Xa — Pa, a/ Eapp\, a 



^A 



(A23) 

(A24) 
(A25) 

(A26) 



with and as defined in Eq. (A4). Analogous to the derivation given by JacksonO for a point dipole, we can 
derive the electric field of a line-dipole at p = OtB 

E(p) = (M(p) - 2nSip)l) ■ p 



(A27) 



where the convention in Eq. (A27) is that the field within the line-dipole at p = is given by the term —2t:S{p)p 
and the Cauchy principal value of the integral should be taken in integrals acro ss the singularity at p = 0. The 
macroscopic field is given by the average over the microscopic field of Eq. (A23) for general p including the chain at 
the origin. Note that since, by symmetry, pA,b ~ —PB,b, the b components do not contribute to the macroscopic field. 
Also by symmetry, we have p^^a = PB,a = Pa- We then have for the macroscopic field Eaifi) ■ 



Ea{p) = -EappKa + 1™ -^P / (P Maa{p - Pk) COs{kUk)pa 

fc^O ZA Junit cell 



27r 

1 



■Pa 



iSappU+hm-T' 



27r 

(fpMaaip) COs{ku)pa ~ -^Pa 



appl,a 



A ■ 



(A28) 

(A29) 
(A30) 



Substituting this result in Eq. ([AT] ) and using the fact that Pa — Pa/A, we find Eq. (|31|): 



1 - — Y ■ 

-L A Xa 



(A31) 



A similar result is obtained for et- 
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3. Evaluation of and 



ab = Mam = = Mrs^ab = Mbm = and 



From the symmetry of Eqs. ( |A22[ ), ( [A24| ) and ( |A25[ ), we see that M 
M^Aaa — ~M^Abb ^^'^ kLs aa — ~M.B bb- This leaves US with Only one element of each matrix to be determined 
Considering first, we split the summation of Eq. ( A24) into two parts. For pi < R (with R large) we perform the 
summation explicitly (taking fc = 0), while for pi > R we replace the summation by an integral: 



M 



■A.,aa 



PieA pl/O pi<R 



k^o 2A 



oc p2Tr 

pdp d(pM_^ 

R Jo 



Xp) COs(fcpCOS (/)), 



(A32) 



which is exact in the limit R —> oo. The sum is evaluated numerically; its values is —0.009677 Aq in the limit R —> oo. 
The integral becomes —tt/{2A) after first taking the limit fc ^ and then the limit i? — > cxj. We can calculate M g 

in a similar way. The sum yields 0.012035 ^ and the integral becomes again — tt/ {2A). Therefore, and M r are: 



M. 



-0.030068 an ^ 



0.030068 an 



^j- _ , -0.008357 



0.008357 an 



(A33) 
(A34) 
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Hntcr (a.U.) 


intra 


8.0 


intra+intcr 
10.0 


12.0 


experiment 




3.59 


2.32 


2.49 


2.69 






1.85 


0.64 


0.76 


0.86 




Eo 


1.74 


1.68 


1.73 


1.73 


1.8 


^Bu — > ^Bu 


0.51 


0.34 


0.39 


0.45 


0.45 


'B, ^ ^Ag 


0.89 


0.45 


0.53 


0.58 


0.55 


TABLE L Quasi-particle [Eg] and optical (-Eo) gaps and binding energies {Eh), for the cases 'intra', using intrachain screening 
only (isolated chain), and 'intra-|- inter', using both intra-and interchain screening (bulk) for three different values of the cut-off 
distance rintcr (see Eq. (j22|)). Exciton transition energies are also listed. Experimental data from Ref. All data in eV. 



Polarizabilities Dielectric constants 



Direction j 




Crystal Axis 7 


e-f{uj — 0) 


1 


60.4 


c 


10.8 


2 


16.3 


h 


3.3 


3 


8.1 


a 


2.6 



TABLE H. The zero-frequency polarizabilities of the single chain per unit length (in ag) and the dielectric constants of the 
bulk along the principal axes of the chain and the crystal. 
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FIG. 1. Monomer and crystal structure of polythiophene. Top panel: A single chain of polythiophene. Black atoms are 

carbon, gray sulfur, white hydrogen. Bottom panel: crystal structure (from Ref. \29\) as seen perpendicular to the chain direction 

(a = 10.5 a.u., b = 14.5 a.u., c = 14.8 a.u., a = 31.2°). 
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FIG. 2. The quasi-particle band structure of both an isolated chain of PT (full lines, left) and bulk PT (full lines, right) 
compared to the DFT-LDA band structure (dashed lines in both pictures). 
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FIG. 3. The hole probabihty along the chain with the electron fixed at 1 a.u. above the inversion center for the two lowest 
singlet excitons, ^B^ (above axis) and ^Ag (below axis), for both the isolated chain (top) and the bulk situation (bottom). 
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